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I discuss the digital disorder introduced by Rosu and Canessa [Phys. Rev. E 47, R3818 (1993)] in the Davydov model of 
energy diffusion along a-helix protein chains. The digitally disordered Davydov beta kinks display self-organized features, i.e., 
power law correlations both in time and space, that may be attributed to incipient dynamic structural changes of the protein 
chain as a consequence of coarse-graining the fluctuations due to microscopic degrees of freedom. In this paper, I provide a 
simple semiconductor model for the flicker noise and also comment on the multifractality that one may associate with the 
protein chains by means of digital disorder. 
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I. INTRODUCTION 
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In a previous paper Q, Canessa and the present au- 
thor, stimulated by the findings of Voss concerning 1/f 
noise in long samples of DNA sequences from GenBank 
[2|, obtained that type of noise in the Davydov's model 
[3| through a simple digital dynamics of a-helix chains. 
I recall that digital dynamics/disorder, i.e., dynamics in 
digital time can also lead to other types of complex be- 
havior, as for example to broken symmetries, oscillations, 
and chaos Q . The aim of the paper is to present possible 
interpretations of the chain dynamics as affected by the 
simple digital dynamics of Q and discuss in a heuristic 
manner the emerging physical picture. 

The outline of the paper is as follows. After reviewing 
the Davydov model and solutions and shortly discussing 
the digital noise of the so-called /3-kink in the next sec- 
tion, I provide in section 3 a simple physical picture for 
the flicker noise in terms of a partial site-trapping of the 
kink. In section 4, I discuss heuristically the kink multi- 
fractality, and I end up with some concluding remarks. 



II. DAVYDOV MODEL 

The a-helix is the most common secondary structure 
of proteins entailing three spines in the longitudinal z- 
direction that we consider infinite in extent and having 

the peptide sequence ( H — N — C = 0)n, n = oo, 

where the dashed lines represent hydrogen bonds. There 
are side radicals and their order is characteristic to each 
protein. The 3 spines are weaved together in a sort of 
incommensurate quasi-one dimensional structure, but 



here we shall consider only the simple one-spine chain, 
containing the peptide groups as molecular units. 

In the Davydov model the main assumption is that a 
substantial part of the chemical energy ( e = 0.422 eV) 
released in the hydrolysis of adenosine triphosphate 
(ATP) turns into vibrational energy (eq — 0.205 eV) 
of the self-trapped amide-I (C=0 stretching) mode of 
the peptide unit. The amplitude of the amide mode 
is self-trapped in the form of a sech- envelope soliton 
whenever there is a balance between the dipolar nearest- 
neighbour interaction and an admittedly quite strong 
nonlinear amide-phonon interaction. The energy transfer 
in the standard Davydov model is through the nonradia- 
tive resonant dipole- dipole interaction and a rather am- 
biguous vibrational-acoustic coherent state, the so-called 
Di ansatz. It is a mixture of quantum and classical 
Hamiltonian methods that has been deeply scrutinized 
in the literature Davydov Hamiltonian can be writ- 
ten in the form 

Hd = Hc=o + Hph + Hint (1) 
Hc=o = 'oBlBn - J(i?iS„+i + (2) 

n 

Hp^,^Y.^n,ib% + ^) (3) 

^^T. x(9)e^'^"^Sii3„(6, -t- 61,) (4) 

The capital and small operators are vibrational and 
phonon ones, respectively. Brown has shown in a clear 
way that Davydov Hamiltonian is a particular case of the 
general Frohlich Hamiltonian of polaronic systems Q 

Hp = Y J„,„at^a„ + ^ fiujqb\bq + ^ hujqixlbl + Xnbq)alan 



qn 



(5) 
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with Di states satisfying the Schrodinger equation of the 
Frohlich Hamiltonian in the hmit Jmn = and another 
D2 ansatz vaHd for Schrodinger evolution in the hmit 



Xr 



0. 



The continuous Hmit, nonhnear Schrodinger (NLS) 
subsonic sohton solutions of the energy transport com- 
ing out from the Davydov model are 
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where ^ = x — Vst is the moving frame coordinate, J 
is the hopping (dipolc-dipole) constant, x is the nonlin- 
ear dipolc-phonon coupling parameter, 7s = l/Vl — 
(s = Vs/va) is the soliton 'relativistic' factor, w is the 
elasticity constant of the chain, ji = y^^l/Jw, Q — 
MRx^ll/2wh^, Es ^ eo-2J + h^v^jAJR^-J^i^/3, and 
eo ~ 0.205el^, the quantum energy of the amide dipole 
oscillator. The first soliton is the vibrational soliton in 
which one may remark the soliton energy Es self-trapped 
by the carrier wave. The second solution is related to the 
local deformation produced by the vibrational soliton in 
the lattice. The /3-kink is a domain- wall configuration of 
the displacements of the peptide groups from their equi- 
librium positions. It has an enhanced stability of topo- 
logical origin since all the peptide groups from the right 
side of the kink > ) are in nondisplaced positions 
whereas all the peptide groups at the left (i^ < 0) arc 

displaced by the same amount f3o — ■ In order to de- 
stroy the domain wall configuration, one should first turn 
the left peptide groups to their initial position. We think 
of the Davydov /3-kink as an interphase boundary for the 
non-equilibrium transition from the Davydov dynamical 
regime of the polypeptide chain to a dynamic 'ferroelec- 
tric' phase of the chain. In the literature on ferroelectric- 
ity it is common to consider the interfacial boundary as 
the kink solution of a time-dependent Ginzburg-Landau 
(GL) equation Our interpretation is based on the 
fact that for low subsonic regime (s^ ^ 1) the Davy- 
dov kink is just the complement of the GL kink (i.e., 
Ku oc (1 — Kql)). Since we are in a non-equilibrium sit- 
uation the more precise terminology for these kinks is dy- 
namic interphase boundaries or interfacial patterns. One 
would like to study their morphology during the growth. 
This is a difficult task since we are in a more complicated 
case as compared to the simple solid on solid model, or 
the kinetic Ising model, where the width of the inter- 
face is given in terms of the nearest neighbor exchange 
interaction ^ , corresponding to the J parameter in the 
Davydov model. The width of the /3-kink is determined 
by two parameters, namely the J of nearest neighbors 
and the x parameter of the nonlinear interaction, con- 
troling the interfacial morphology. 



Perhaps, one should notice the formal analogy between 
the form of the Davydov kink and the Glauber transition 
rate in one-dimensional spin chains. The message of this 
analogy is that the kink is just a step structural function 
which is required by a Hamiltonian evolution and by a 
detailed balance condition in the spatial coordinate. 

The digital disorder introduced in Q| is due to small 
random displacements of the instantaneous kink position, 
which a priori may appear to be equivalent to small fluc- 
tuations of the kink velocity. The procedure is as follows. 
The centre of mass of the kink in the moving frame is 
taken as the spatial origin and small random excursions 
Aa; are allowed between ±1. The time correlations of 
such fluctuations are estimated by means of the dimen- 
sionless noise power spectrum, which usually characterize 
random processes of the quantity chosen as probability 
density, in our case the (3- kink function. Thus the tem- 
poral correlations at the scale < t <r«l// are 
measured by B 



T Jo 



(9) 



When calculated with a standard fast Fourier transform 
(FFT) algorithm, S/sif) displayed clearly 1/f noise [|. 
Moreover, the fact that the hyperbolic tangent is the 
solution of the GL equation with real coefficients leads 
to power law correlations in the space domain through 
a multifractal formalism. But power law correlations 
both in time and space are necessary requisites of self- 
organized criticality (SOC) |10|. The digital noise has 
its origin in the removed microscopic degrees of free- 
dom, for instance the type of lattice, when one is passing 
to the continuous NLS limit. In this sense, the digital 
noise is an internal one in the terminology of the SOC 
literature while the disorders coming from varying 
the parameters in the Davydov model, which were stud- 
ied some time ago by Forner belong to the class 
of external noises. Dfaz-Guilera showed that the 
two types of noises can be considered as different limits 
of an Ornstein-Uhlenbeck stochastic process. Actually, 
the digital noise is very close to a proposal of Cruzeiro- 
Hansson [|l3|, who has argued that the thermal desta- 
bilization of the Davydov soliton at low temperatures 
starts through a cascade of disordered states in the dis- 
placements, of intermediate energy between the soliton 
and the exciton states. At higher (biological) tempera- 
tures, the disorder in the displacements is transmitted 
to the vibrational soliton which in this way is destroyed, 
although the states populated in the decay preserve the 
localization of the amide I excitation [Q. As a mat- 
ter of fact, experimental evidence provided by optical 
calorimetry shows that phonons and vibrations heat up 
at different rates [[151. 
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III. MODEL FOR THE FLICKER NOISE OF THE 
/3-KINK 

A simple interpretation of 1/f noise may be given in 
terms of a partial site-trapping of the /3-kink. This trap- 
ping modulates the mobihty of the domain wall and is 
correlated with the carrier- wave self-trapping. The toy 
model is adapted from semiconductor physics ||l^. We 
shall consider the a-chain to be at a slightly non-zero 
temperature, say T < IIK |l^], where the Davydov sys- 
tem is still not washed out by thermal effects. Let us sup- 
pose a sort of site- trapping changing the kink mobility 
by rectangular pulses, whose duration r is the trapping 
time. Furthermore, all pulses can be considered statisti- 
cally independent and thus their distribution is Poisson 
p{t) = r exp(— r • t), where t is the time interval between 
succesive pulses and r is the average number of pulses per 
unit time. Next we shall consider the distribution func- 
tion of the trapping times r, due to a single site-trap to be 
of the form f;(r) = exp(— t/tq), where tq is the average 
trapping time. The last assumption is that different site- 
traps have different tq's, which are related to the trap- 
ping energy E in the thermal way, tq = Kexp{—E/kT), 
where K is a constant, and E is distributed among the 
trap ensemble according to a square distribution spec- 
trum between two limiting energies Ei and i?2. From 
the thermal assumption one obtains 



1 



1 



log(r2/ri) To 



(10) 



valid for ti < tq < T2, where ti = if exp(— _Bi/fcr) 
and T2 — exjp(~E2/kT). The power spectrum of a time 
series of rectangular pulses with the same duration t and 
Poisson distributed, was obtained long ago by Gisolf jl^] 
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By combining ^(t) with Eq.(ll), Van der Ziel ||l^ ob- 
tained a shot-noise power spectrum 



Si', — const ■ 



To 



(12) 



One gets the 1/f noise by averaging the shot spectrum 
with the square distribution of the trapping times Eq. (llC 
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kT 1 



flicker 
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Of course, the simple flicker mechanism we have used 
here was more for illustrative purposes, and further work 
is required to clarify the nature of such noise in protein 
chains. For example, the flicker noise can be explained 
alternatively by means of a scaled Langevin equation as 



introduced by Kogama and Hara For that, it is 

sufficient to consider the same scaling parameter both in 
the Langevin variable and in the attached noise as was 
shown in their paper. 



IV. MULTIFRACTALITY OF /3- KINKS 

The implicit occurence of the GL kink in the Davydov 
(3 displacements is essential for discussing the multifrac- 
tality issue. Previously, Brax |2^] showed the equivalence 
of the GL equation with real coefficients under random 
initial conditions and the linear heat equation with Gaus- 
sian random potential and made a multifractal analysis of 
the problem of direct relevance for our study. As is well- 
known the GL equation, which is a cubic reaction- di- 
fussion equation, describes phenomenologically the evo- 
lution of the order parameter in superconductive phase 
transitions, and models also spatial and time fluctuations 
of systems near Hopf bifurcations. In 1972, Scalapino, 
Sears and Ferrell ||2^ studied in detail the statistical me- 
chanics of one-dimensional GL fields. They remarked 
that such fields can describe the dynamical behaviour of 
nearly-ordered systems which are not undergoing sharp 
phase transitions, and conjectured that the real-field case 
may have application in some organic chain systems. Our 
approach is just a further confirmation of that idea. As a 
matter of fact, in the D'Alambert variable {S^ = x — Vgt), 
the GL kink Kq l is the solution of a GL equation of the 
type 



GL 



dt 



^d^KcL/dx^ +PiKgl 



(14) 



with real pi and coefficients. Such GL equations are 
typical for the structural phase-transitions H] in equilib- 
rium situations, but they can be used also in nonequilib- 
rium/drivcn systems. Moreover, in Eq.(8) one can recog- 
nize the complement of the GL kink (i.e., the Davydov (3- 
kink) playing the role of the probability density distribu- 
tion . Following Brax, one can develop a multifractal 
formalism if the kink probabilistic distribution function 
is identified with the r- function of the multifractal for- 
malism 



G{^2)-G{^i) 



(15) 
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1 recall that in the multifractal formalism the function 
Tq is the cumulant generating function. The derivative of 
this function with respect to q is denoted by a and via the 
Lagrange multipliers procedure of statistical thermody- 
namics, one obtains the function f{a) which is the den- 
sity of a measure, and could be interpreted as a fractal 
dimension when it is positive, and related to instabili- 
ties for negative values. The equations a — Sr/Sq and 
/(a) = qa — T represent the basis of the multifractal for- 
malism pJ. Formally, q is the inverse temperature, r 
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is the Gibbs free energy, and / is the entropy. One can 
plot the second derivative of the r-function (the 'specific 
heat') and find out intervals of the 'temperature' vari- 
able, within which the plot clearly displays features of a 
phase transition, that is a peak in the 'specific heat' at a 
certain value of the 'temperature' variable 

V. CONCLUDING REMARKS 

The digitally disordered Davydov beta kinks can dis- 
play features characteristic to self-organized criticality. 
Here, I considered the digital noise of Davydov kinks as 
a reflection of hidden microscopic degrees of freedom of 
the protein chain. It is an internal noise in the SOC ter- 
minology, representing the dynamics of a self-organized 
state of the chain, therefore not destroying the soliton 
solution all over the scaling region. 

I ascribed the 1/f noise to carrier trapping and so to 
some sort of short-distance thermal activation physics, 
rather than to coherence effects. This would be more 
in the standard spirit of solid state physics. In general 
1/f mesoscopic noise can have many interpretations and 
it would be helpful to disentangle its real origin. At the 
level of vibrational degrees of freedom it will be of interest 
to see if the vibrational soliton possesses self-organized 
features as well, and to investigate the connection with 
the kink digital disorder. 
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